On the matrix elements of AB = operators 
in the heavy meson decay widths 
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Q ' Abstract 



We determine the chiral corrections to the matrix elements of the AB = four- 
quark operators which are relevant to the studies of the ratios of lifetimes of heavy- 
t:;j- ■ light mesons as well as to the power corrections to the inclusive semileptonic heavy- 

C^ , to-light decays. The chiral logarithmic corrections computed here can be combined 

with the forthcoming estimates of the corresponding matrix elements on the lattice 
to provide the reliable physics result of the well known bag-parameters Bi^2 and £1^2- 
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1. Phenomenological introduction: The matrix elements of AB = operators enter 
several phenomenological studies of which the most important ones are the analyses of the 
spectra of inclusive semileptonic decays of heavy mesons and the lifetime ratios of heavy- 
light mesons. 

1.1. Power correction to r{B -^ XuCu): Controlling the power corrections in the 
spectra of inclusive semileptonic heavy to light decays has been -and still is- an important 
obstacle when aiming at the reliable extraction of the corresponding CKM parameters [1]. 
This is particularly important in the case of \Vub\- In ref. |2] it has been shown that the 
l/?Ti^-corrections involve the matrix elements of dimension-6 four-quark operators of the 
flavour structure AB = 0. More specifically 



where the matrix elements are conveniently expressed in terms of bag parameters, -81,2, 
as [3] 



{B\0^s-p\B) ^ {B\bil - ^,)uu{l - ^,)b\B) = flmlB, 



(2) 



with fs being the S-meson decay constant. Therefore what is actually needed in eq. ([T]) 
is the difference of bag parameters, Bi — i?2- The early estimates of B12 in the framework 
of QCD sum rules in the static heavy quark limit were reported in ref. [1], and recently 
extended to the full QCD case [5j. They were also computed on the lattice. In ref. [6] 
the authors found that in the static heavy quark limit Bi — B2 is zero, while the lattice 
study with the propagating heavy quark indicated that Bi — B2 can be quite different from 
zero [7]. We are now in the era of ever better unquenched lattice QCD studies and a new 
lattice computation of the bag parameters of the AB = operators is clearly desired. In 
recent years it became evident that the control over the chiral extrapolation is essential in 
order to reduce the systematic uncertainties in the results of the lattice QCD studies. In 
this paper we provide the chiral corrections associated with the bag parameters computed 
in the static heavy quark limit. 

1.2. Bi — B2 from D^s: Before we turn to the question of chiral corrections, let us see if 
we can get some information about the size of Bi — B2 from the available information on 
D-decays. The expression for the inclusive semileptonic decay width up to and including 
the terms oc 1/ml, and neglecting the small contribution oc \Vcd\ can be written as 
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where the phase space factor 

Iq(z) = {l-z'^){l-8z + z^)-12zHogz, h{z) 

and z = ml/rn^^, while the numerical parameterisation of the a^-correction to the partonic 
decay width r][z) reads [5], 



^r 



(4) 



7]{Z) 



2a, 



TX 



^''. 



z^'-f + l 



The equation of motion allows us to write 

^ (ciiDA^c + c^a.G(^ + C(l/m 
\mj. V 2 / 



cc = ape 



cl -I 



which in the standard notation ^ 



l^l 



27n 



D 



-{D\c{iD^fc\D) 



f^G 



-^(D\c^a.Bc\D) 



(5) 



(6) 



(7) 



^Also standard is the notation in terms of Ai 2, the parameters measuring the kinetic and chromomag- 
netic energy of the heavy quark inside a heavy-Ught system. The relation to yUg ^ is: Ai = — yU^, and 
A2 = Mg/3- 



can be written as 



1 (flNfl) = l-4-#. (8) 
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Finally eq. ([3]) becomes 



TiD^Xeu) = ^|K.|%(;^M.)|l + ^^[/^^-(l-4/i(.)//o(.))/.a| 

-^^\V^s?r^{z)flmn{B,-B,) . (9) 

Clearly the heavy quark expansion applied to the decay of charmed mesons is expected 
to converge much slower than in the case of 5-mesons. It is however interesting to use 
the available information on charmed modes to bound the Bi — B2 value. Concerning the 
parameters appearing in the first line of eq. ([9]) we can use jj^q = (3/4) [m|,. — m|,] = 
0.41 GeV^, while the value of /i^ is still somewhat vague. Recent experimental fits to the 
moments of the semileptonic b —>■ c decay spectrum P, [10, [H] quote rric ~ 1-1(1) GeV 
and /i^ ~ 0.5(1) GeV^ in the so-called kinetic scheme [12] and at yU = 1 GeV. When 
converted to the MS scheme, the charm quark mass is mcirric) = 1-2(1) GeV, consistent 
with the estimates based on the lattice QCD simulations mc{mc) = 1.32(3) GeV, and 
1.30(3) GeV [13], as well as with the recent QCD sum rule study mdmc) = 1.29(1) GeV [1] 
If, in addition, we take fn = 208(4) MeV [IS], r^± = 1.040(7) ps, r^o = 0.410(15) ps 
and the recently measured semileptonic branching fractions [T7]: ^ 

B{D^ -^ Xeu) = (16.13 ± 0.20 ± 0.33)% , B{D" -^ Xeu) = (6.46 ± 0.17 ± 0.13)% , (10) 

then we get that the terms including 1/m^-corrections saturate the experimental value for 
the branching ratio to about 85%. If we now assume the 1/ml term fully saturates the rate, 
we get B2{mc) — Bi^nic) ~ 0.04. When evolving those bag parameters to the m;,-scale [5|[TU]. 
the difference is further increased by about 30%, leading us to B2{mb) — Bi{mb) ~ 0.05. 
After plugging that number in eq. ([1]), and by taking /b = 0.2 GeV, we get 



B{B ^ X„ez/)i/„,3 < 0.04 x IK^I', (11) 

which is {comfortably) a very small number. Of course this exercise is only a speculation, 
while for the reliable estimate of B{B — »■ Xuez/)]^/^3 a direct non-perturbative method 
should be employed to compute the matrix elements |2]). 

1.3. i?-meson lifetimes: When studying the 5-meson lifetimes, due to the nonleptonic 
decay modes two more operators enter the game 

{B\T^^^\B) ^ (5,16^7,(1 - l5)qqY^'{l - l^)b\B,) = /|m|ei , 

{B\T^^^\B) ^ {B,\b^{l-^,)qq^{l-^,)b\B,) = f^mle^ , (12) 



^Very recently the other charm factory (BES) presented similar results but with the final muon instead 
of electron [18]. Their values are fully consistent with those given in eq. pO)) . measured by CLEOc, but 
the error bars are an order of magnitude larger. 



where A"^ are the Gell-Mann matrices. The bag parameters in eqs. (l2][T2l) were first in- 
troduced in this form in ref. [3] in which the authors studied the lifetime difference of 
hadrons containing one valence fe-quark. Voloshin, however, reahsed in ref. [2] that the 
light flavour content of the operators can be different from the light valence quark of the B- 
meson. Those are the (in)famous "eye-contractions" which are extremely difficult to study 
non-perturbatively. In addition, an extra penguin operator contribution was singled out in 
ref. [19]. The contribution of the non- valence and the penguin operators are expected to 
be negligible in the case of the meson lifetime ratios due to the light flavour symmetry of 
the spectator quark. That argument however does not apply to the power correction to the 
semileptonic decays ([3]). In terms of bag parameters and neglecting the eye-contractions as 
well as the penguin operator contributions, but including the NLO QCD-corrections to the 
Wilson coefficients computed in refs. [19l [20] , the master formulas for the lifetime ratios of 
5-mesons read 



r{B,) 

rjBs] 
r{B,) 



1 + 0.07(2) X Bi + 0.011(3) X 52 - 0.7(2) xei + 0.18(5) x €2 , 

1 + 0.007(2) X [Bl - Bi] - 0.009(2) x [B'^ - O.9B2] 

+0.15(4) X [el - l.l£i] - 0.18(5) x [e^ - 0.9ei] , (13) 



where the superscript "s" has been used to distinguish the bag parameters for the case of 
valence strange quark. In this case it is even more important to have a good handle on the 
£1,2 parameters whose impact is enhanced by the size of the Wilson coefficients (numerical 
values of which are displayed above). 

2. Chiral corrections : One of the main problems in relating the bag parameters B12 
and £i_2 computed on the lattice to the physical bag parameters is the necessity to perform 
the chiral extrapolation of matrix elements computed with the light quark masses directly 
accessible on the lattice (1 > mq/ms > 1/4) down to the physical limit {mq/nis ~ 1/25). ^ 
The expressions derived in chiral perturbation theory provide an important guidance in 
that respect. The chiral corrections to the matrix elements of the whole basis of four-quark 
AB = 2 operators were recently computed in refs. pTl [22] . We showed in ref. [22] that the 
validity of the formulas derived in heavy meson chiral perturbation theory (HMChPT) may 
be questionable for the quarks not lighter than about the third of the strange quark mass, 
because of the nearness of the scalar heavy-light mesons (or more precisely, of the heavy- 
light (|)"'"-doublet). In other words, unless one wants to deal with a very large number of 
low energy constants, the adequate HMChPT expressions are only those with Nf = 2 light 
quark flavours (i.e. with the pion loops only). In this paper we do not return to that issue. 
Instead we focus on the chiral corrections to the bag parameters of the AB = operators 
introduced above. 

2.1. Pramework: As in ref. [22], in the present paper we work in the static heavy quark 
limit and use the HMChPT lagrangian already described in detail in ref. [22]- We chose a 



^Here and in the following niq = iriu = md- 



basis of operators 

Ov-A = &7/.(l -75)gLgL7''(l-75)&, 

Os-p = 6(1 -75)gL 9^(1 + 75)6, 
Ts-P = 6(1 -75)tVgLi^(l + 75)6 • 

The heavy quark spin (S) and the chiral symmetry {Ul^Ur) transformations act on the 
heavy and hght quark respectively hke, b -^ Sb (i.e., ■yob = b), and qL,R -^ UL,RqL,R- 
Since the colour structure [short distance) does not influence the chiral logarithms [long 
distance) [23], we need to consider only two of the above operators which we choose to 
be Ov-A and Os-p- In HMChPT we need the bosonised forms of these operators, which 

are built up from the heavy-light (i)~-doublet flelds, Hq{v) = -^ \P'^J^v)y^ — Pq(w)75] , 
and the pseudo-Goldstone flelds, S = exp(2i0//), where is the usual matrix of pseudo- 
Goldstone bosons. Under the heavy quark and chiral symmetry the fleld Hq{v) transforms 
as Hq -^ SHqiU]^,^, while S transforms as S ^ Ul^Uj^. The standard procedure then 
consists in introducing ^ = yS = exp{i(l)/f), which transforms as ^ — * Ul^U^ = U^Uj^. 
By a simple chiral and heavy quark spurion analysis of the operators Oy-A^ Os~p we then 
obtain their most general bosonized form within HMChPT, namely 

0^_^ = Y. ^i^Tr [m),l,{l - 1,)X] Tr [^7^1 - l,)m\] 

X 

+ Y, ^i^Tr \Hqn,{l - 75)X] Tr [X7'^(l - y,)Hq,] + c.t. , 

X,q' 

OLp = $^r2xTr[(e:H^),(l-75)X]Tr[X(l + 75)(i/e^),] 

X 

+ Y ^2xTr [Hq,{l - 75)X] Tr [X(l + y,)Hq>] + c.t. , (14) 

X,q' 

where "c.t." stands for counterterms, and X G {1,75,7^,71/75, cr^p}. The contraction of 
Lorentz indices and HQET parity conservation requires the same X to appear in both 
traces in the products. Any insertion of -^ can be absorbed via equation of motion, '^H{v) = 
H{v), while any nonfactorizable contribution with a single trace over Dirac matrices can 
be reduced to the form written in (TT4l) by using the 4x4 matrix identity 

4 Tr{AB) = TT{A)Tr{B) + Tr(75A)Tr(75fi) + Tr(A7^)Tr(7^5) 

+Tr(A7^75)Tr(757^B) + ^TT{Aa,,)TT{a>'''B). (15) 

Note that these matrix identities also ensure invariance of the bosonized operators under 
Fierz transformations, as can be readily checked by bosonizing the Fierz transformed quark 
level operators, which involve formally different spurion flelds. 

An important comment is also that the second lines in the operators written in eq. (IHl) 
stand for the corresponding "eye-contractions" . The sum over q' runs over all light quark 



flavours. Once saturated by the external heavy-light meson states of the light flavour q only 
the eye-contraction with q = q' will contribute. However at one loop in HMChPT these 
contractions will not produce any chiral logarithmic correction to the matrix elements of 
the above operators. Their effect can show up at two or more loops though. 

After calculating the above traces, and by retaining the pseudo-Goldstone fields up 
to quadratic order, we obtain: 
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where, for simplicity, we do not display the counterterms and we used 
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To 



^1 + T175 — 4(ri^„ + n^^^s) — 12ri^^p, 

T"2 — ^275 — ^27^ + ^27^75 + ^2 — ^"275 — ^27„ + ^"27^ 75, 
^27^ — ^27^75 + T2^^ — T2^^^5 . 



(17) 



Similarly S[*2 stand for the combinations of 5ix,2X couplings appearing in eq. 0141) . 
2.2. Chiral loop corrections The computation of the chiral loop corrections to our 
operators is by now standard. It involves 6 diagrams which are shown in fig. [H Four 
graphs are factorizable and two are not. Of factorizable diagrams we have the self energy 
contributions [(a) and (b) in fig. [T] which give rise to the wave function renormalisation 
corrections {SZg) which can be found in e.g. [22], and two tadpole graphs which represent 
the loop corrections to the weak currents composing the four-quark AB = operator 
[(c) and [d) in fig. [1]. Finally there are two nonfactorizable graphs: tadpole [(e)], and 
"sunset" [(/)]. An alerted reader may notice the absence of the mixed terms, i.e. the ones 
involving an exchange of a pseudo-Goldstone boson between the weak operator and the 
HMChPT interaction Lagrangian. Those contributions drop out due to heavy vector meson 
transversality {v ■ Sp* = 0). Compared to the situation we encountered in the computation 
of the chiral loop corrections to the matrix elements of AB = 2 operators [22], in the 
present situation the sum of diagrams (c), [d] and (e) vanishes. 
The resulting expressions read 
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Figure 1: The graphs giving the non-zero contribution to the NLO chiral corrections to the matrix 
elements of the AB = operators discussed in this paper. The double lines correspond to the 
heavy-light mesons, and the dashed ones to the pseudo-Goldstone bosons. AB = operators are 
denoted by '1X1X1", while the strong vertices coming h'om the HMChPT lagrangian are denoted by 
the full dots, "•". 



where JqK) = K/47r)2 log K//i'), Y2 = {B*/B2f''^ X, 
and f are the SU{N) generator matrices. Summation over "i 
understood. 



= (f,/i?,)T- ^ (1 - 6,/n), 
in the above expressions is 



3. Results : On the basis of the expressions derived in the previous section we now discuss 
our results. We will first give the explicit formulas for the chiral corrections that might 
be particularly useful to the lattice practitioners. We will then make some important 
assumptions which will allow us to infer a few phenomenological implications. 
3.1. Message relevant to the lattice QCD studies: In our previous paper we showed 
that due to the nearness of the (|) -doublet of the heavy-light mesons, only the pion loop 
contributions are a safe prediction of this (HMChPT) approach which then can be used 
to guide the chiral extrapolations of the heavy-light meson quantities computed on the 
lattice. This was shown to be the case for the decay constants, Standard Model and SUSY 
bag parameters parameterising the matrix elements of the AB = 2 operators [22] , pionic 
couplings g, g and h [25], as well as for the Isgur-Wise functions [21]. The same holds true 
in this case. Therefore the relevant expressions to be used in the lattice extrapolations in 
the light quark mass are those derived in S'?7(2)-theory and they read: 
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or by recalling that 
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for the bag parameters we have 
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where g'^ can be computed separately on the lattice as in ref . [26] , and the parameters of the 
fit are Bj^^^ and the counterterms &i,2(yu)- It is worth emphasizing that the /i-dependence 
in the chiral logarithms cancels against the one in the low energy constants. The situation 
with the chiral corrections to the matrix element Os-p is similar to what we discussed in 
ref. [22] where for the non-Standard Model AB = 2 operators the new low energy constant 
"F" appeared. Its value is likely to be very close to unity as it represents the following 
ratio 
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and it can be relatively easily evaluated on the lattice. "^ Finally, let us stress once again 
that thanks to the identity 
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the chiral logarithms to the bag parameters £1^2 are of the same as those in B12 parameters 
but their low energy constants are of course different. To be fully explicit: 
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3.2. Back to phenomenology: In the early phenomenological applications the formulas 
derived in HMChPT were used to estimate the size of the hadronic quantities by using 
the theory with Nf = 3 light flavour and by neglecting the counterterms (or, at best, 
estimating them by means of some quark model). Nowadays we also know that the (i) - 

states should be included if one is to use HMChPT with Nf = 3. In what follows, the (|) - 
contributions will be neglected too, which is an extra assumption. To get the difference 
B2—B1 we will proceed along these lines and impose Bj^2'^ = 1, like in the vacuum saturation 
approximation, and neglect the counterterms, to obtain 
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(27) 



This is as far as one can get at this stage, since there is no information available concerning 
the size of Y2. We reiterate that it can be computed on the lattice as indicated in eq. (122|) . 
Note in passing that if we use B2 — B1 = 0.05 as inferred in introduction from the D-decays, 
we would obtain Y2 ~ 0.8. Similarly, for the bag parameters entering eq. (1131) we have 
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(29) 



"'Notice that the complexity related to the matching of the operator computed on the lattice to its 
counterpart renormalised in the continuum renormalisation scheme completely cancels in that ratio. 



9 



which, together with the assumption that ej^2'^ = £^2^ = brings eq. (TT3|) to 



r{B,: 



1.077 + 0.002 Fa (1-076 ± O.OOS)^''" , 

1.001 - 0.002 Fa (0.950 ± O.Oig)^''^ . (30) 



where in the parentheses we also give the experimental values [TB]. We see that in spite of the 
assumptions the current experimental information does not allow to constrain appreciably 
the value of the coupling 1^2- 

4. Summary: In this paper we presented the result of our calculation of the chiral cor- 



rections to the matrix elements of four-quark AB = operators that are relevant to the 
phenomenology of the lifetime ratios of the heavy-light mesons and to the inclusive semilep- 
tonic decay spectra [27]. The calculation of the chiral corrections can be combined with the 
lattice calculations of the AB = matrix elements to either extrapolate the lattice data 
towards the physical light quark masses, and/or to fix the counterterm coefficients 61,2 (/^) 
in eq. fl2Tl) . the couplings Y2 fl22l) and Y2 f l25l) and the tree level bag parameters (in terms 
of chiral expansion) . 
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